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ABSTRACT 

Motivated by upcoming data from astrometric and spectroscopic surveys of the Galaxy, 'we explore the 
chemical abundance properties and phase-space distributions in hierarchically-formed stellar halo simulations 
set in a ACDM Universe. Our sample of Milky-Way type stellar halo simulations result in average metallicities 
that range from [Ee/H] ~ -1.3 to -0.9, 'with the most metal poor halos resulting from accretion histories that 
lack destructive mergers 'with massive (metal rich) satellites. Our stellar halo metallicities increase 'with stellar 
halo mass. The slope of the [Ee/H] -M* trend mimics that of the satellite galaxies that were destroyed to build 
the halos, implying that the relation propagates hierarchically. All simulated halos contain a signihcant fraction 
of old stellar populations accreted more than 10 Gyr ago and in a few cases, some intermediate age populations 
exist. In contrast with the Milky Way, many of our simulated stellar halos contain old stellar populations 
which are metal rich, originating in the early accretion of massive satellites (M, ~ IO^Mq). We suggest that 
the (metal rich) stellar halo of M31 falls into this category, while the more metal poor halo of the Milky Way 
is lacking in early massive accretion events. Interestingly, our hierarchically-formed stellar halos often have 
non-negligible metallicity gradients in both [Ee/H] and [a/Ee]. These gradients extend a few tens of kpc, and 
can be as large as 0.5 dex in [Ee/H] and 0.2 dex in [a/Ee], with the most metal poor halo stars typically buried 
within the central ^ 5 kpc of the galaxy. Thus evidence for metallicity gradients alone in the Milky Way stellar 
halo would not preclude its formation via a hierarchical process. Only coupled with phase-space data can 
metallicity information be utilized to test ideas about accreted versus in situ formation. Einally, we hnd that 
chemical abundances can act as a rough substitute for time of accretion of satellite galaxies and, based on this 
hnding, we propose a criterion for identifying tidal streams spatially by selecting stars with [a/Ee] ratios below 
solar. 

Subject headings: galaxies: abundances — galaxies: evolution — cosmology: theory 


1. INTRODUCTION 

The stars in the halo hold important information about the 
formation history of the Galaxy. Their spatial and velocity 
distributions can be used to retrace their dynamic origin, and 
their chemical abundances can constrain the star formation 
histories of their constituents. A major goal over the next 
decade is to obtain kinematic and chemical information for 
a large number of stars in our own Galaxy, thus providing a 
detailed reconstruction of its fo rmation history and a link to 
the un derlying cosmology (e.g. lEreeman & Bland-Hawthornl 
I2002h . It is the goal of this work to provide a fist step to¬ 
wards the modeling needed to exploit these observations to 
their fullest potential. 

The combination of spatial, kinematic and chemical data 
will provide a powerful discriminate for testing different 
Galaxy formation models. Historically, metallicity gradients 
have acted a primary motivator in formulating hypotheses: 
large-scale metallicity gradi ents are associated with a rapid, 
semi- continuous collapse dEggen, Lvnden-Bell & Sandagel 
[T^ ELS); homogeneous metallicity distribut ions are asso¬ 
ciated with chaotic assembly from fragments JSear a cm 
ISearle & ZinnII 19781 SZ). While there are no well-motivated 
models that predict these opposite extremes from hrst prin¬ 
ciples, the poles of debate set by ELS and SZ serve as im¬ 
portant straw-man theories for shaping ideas and testing data. 
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Nevertheless, realistic models will always fall into less ideal¬ 
ized categories. Eor example, our ACDM-based model relies 
entirely on a hierarchical origin for the stellar halo, yet pre¬ 
dicts some non-negligible metallicity gradients (see below). 
Likewise models that allow significant “rapid” in situ star for¬ 
mation necessarily contain some accreted stellar material as 
a result of having a background hierarchical cosmology (e.g. 
iRenda et alJl2005t) . Clearly, only in combination with spatial 
and kinematic data can chemical abundance data be used to 
its full potential. 

Numerous results over the past decades support the va- 
lidity of the hierarc h ical model of s t ructur e formation 
jWhite & ReesI l1978t iRliimenthal et al.l 11984) . In this 
model, the expectation for stellar halos of galaxies like 
the Milky Way is that they form in large part by phase¬ 
mixing of tidal debris from the numerous accreted satel¬ 
lites. Although the dynamical evolution of galaxy halos 
has been extensively modeled, and we have now predic¬ 
tions for the phase-space distribution of b oth dark matter 
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Intuitively, one expects stochastic accretions to leave 
peculiar signatures in both the kinematics and chem¬ 
ical abundances of present day stars. Some obser¬ 
vational evidence in that respect is indeed found in 
our Galaxy. The halo contains stars which stand 
out in both metallicity and velocity JCarnev et alJ 1 19961 
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stellar stream in M31, a merger debris extending more than 
100 kpc to the SE of M31’s disk, stands out from the back- 
ground halo not only in stellar over-densities and k inematics 
[ibata et al.l2001all2004lGuhathakurta et al.l2005bl) . but also 
in metallicity dFerguson et al.l20()^. 

Over the next decade an immense amount of data will 
become available for stars in the Galaxy from both astro- 
metric and spectroscopi c surveys - e.g. from the astromet¬ 
ric satellite GAIA (e.g. iPerrvman et S ] 1200ll ), the “R adial 
Velocity Experiment” (RAVEl (e.g. l■Steinmet7^l2003^ ■ and 
the “Sloan Extension for Galactic Underpinnings and Ex¬ 
ploration”, segue'. Positions, line of sight velocities and 
proper motions will be measured for millions of stars in the 
halo, enabling the full phase-space to be reconstructed. Simi¬ 
larly, high accuracy, wide field measurements of stellar spec¬ 
tra will allow the mapping of the Galaxy in chemical abun¬ 
dances. These studies will not only constrain cosmology and 
galaxy formation (e.g. lBullock, Kravtsov. & Weinbergi2001l) 
but also test ideas a bout first light and reionization (e.g. 
iTumlinson et al.l2004l) With the upcoming observational data, 
the goal of putting together the history of our Galaxy is fi¬ 
nally becoming feasible. It is therefore important to have the¬ 
oretical predictions for the combined kinematic and chemical 
abundance distributions, as well as to have a coherent theoret¬ 
ical framework for interpreting the upcoming results. 

Although future surveys promise to provide detailed maps 
of the chemical and phase-space structure of the Galaxy, we 
may worry that ours is but one among many galaxies of its 
size. How typical is our galaxy? How can it be used to con¬ 
strain general ideas about cosmology and galaxy formation 
if it is atypical in some way? Fortunately, in the context of 
ACDM, we have well defined expectations for variations in 
the formation times and accretion histories of Milky-Way size 
galaxi es and this idea is at the hea rt of our exploration. For ex¬ 
ample ll200.5albh recently have estimated the 

metallicities of the inner halos of several nearby spiral galax¬ 
ies and find that they are quite metal rich compared to the 
stellar halo of the Milky Way. They go on to suggest that 
the Milky Way is not typical for a normal spiral galaxy of 
its luminosity. As discussed below, if this result holds then 
our models provide a straightforward interpretation: that the 
Milky Way has experienced fewer than average major accre¬ 
tion events with LMC-size objects. Indeed, phase space infor¬ 
mation could test ideas of this kind specifically by revealing 
few major accretion events and many more lower-mass accre¬ 
tions. 

The modeling of the chemical evolution of the Galaxy in a 
reali stic cosmological context is just beginning to be explored 
(e.g. lRekki & Chiball200lHRrook et alji2003l 120041) . Recent 
advances in modeling techniques, like the advent of “hybrid” 
methods, allow us to address this problem with greater reso¬ 
lution. Hybrid methods enable the modeling of star formation 
and chemical enrichment with detailed semi-analytical pre- 
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scriptions, and at the same time provide th e detailed dynam¬ 
ics vi a the coupled N-body simulations (e.g. lKaiiffmann et alJ 
Il999t ISpringel et aljEoOll) . In this study we use a previ¬ 
ously developed hybrid method llBulIock & Johnstonll2005t 
[Robertson et al.ll200.^ [Font et alJl200-5l) to explore the phase- 
space distribution of chemical abundances in a series of Milky 
Way-type galaxy halos formed in a ACDM Universe (H,„ = 
0.3, Ha = 0.7, h = 0.7, and ag = 0.9). 

A description of our sample of stellar galaxy models is 
given in Section |2] In Section 0 we present results on the 
spatial distribution of [Fe/H] and [a/Fe] chemical abundances 
in these stellar halos and propose a criterion for detecting cold 
stellar streams based on their [a/Fe] abundances. In Section 
l^we discuss our results and compare them with the available 
observations and in Section|3we conclude. 

2. THE SAMPLE OF MILKY WAY-TYPE STELLAR HALOS 

In a previous study iRullock & .lohnst^ li2005l) have pre¬ 
sented a sample of 11 stellar halos of Milky Way-type galax¬ 
ies assembled hierarchically from accreted satellite galaxies 
in a ACDM cosmology. The physical properties of the simu¬ 
lated halos are summarized in Table 1. The merger histories 
of these halos have been specifically selected not to have high 
mass recent mergers, so as to maximize the pro bability that i t 
will host a disk galaxy like the Milky Way (e.g. lWvsal200l[) . 
Our model stellar halos have similar mass, density profile and 
total luminosity as the Galactic halo. The models also obtain 
roughly the same number of surviving satellite galaxies as in 
the Milky Way and match their physical properties: the stel¬ 
lar mass-circular velocity relation, M*-Vcirc, for Local Group 
dwarfs; and the distribution of surviving satellites in luminos¬ 
ity, central surface brightness and central velocity dispersion. 

In this paper we extend our models to include chemi cal evo¬ 
lution with the prescriptions of iRobertson et ail ll2005l) . These 
include enrichment from both Type II and Type la supernovae 
and feedback provided by supernovae blow-out and winds 
from intermediate mass stars. We follow the chemical evo¬ 
lution of each satellite galaxy accreted onto the main Galaxy 
and trace the build-up of a series of metals, such as Ee or a- 
elements. Each star particle in the simulations has assigned a 
set of [Ee/H] and [a/Fe] abundance^ ratios determined by the 
star formation history of the satellite in which the star orig¬ 
inated. The coupled N-body simulations resolve each Milky 
Way-type stellar halo with a few million particles and allow 
us to explore the phase-space distribution of chemical abun¬ 
dances in unprecedented detail. 

3. RESULTS 

3.1. Global Properties of the Simulated Halos 

Tidal streams in the halo display various degrees of spatial 
coherence, given their different times of accretion and orbital 
properties. Similarly, we expect to see various degrees of co¬ 
herence in the chemical abundance space, given the different 
star formation and chemical enrichment histories of their pro¬ 
genitor satellites. 

Figure ^ shows the present day (f = 0) spatial distribu¬ 
tion of [Fe/H] and [a/Fe] abundance ratios for halos HI - 
H8. The maps span 300 kpc on a side {x and z directions, 
respectively) and 300 kpc in depth (y). Each pixel repre¬ 
sents a color coded chemical abundance, calculated as the 
average [Ee/H] or [a/Ee] of all stars included in a volume 

^ The [a/Fe] ratio is defined as the average of [Mg/Fe] and [0/Fe] ratios, 
and used thereafter in the paper. 
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Halo 

# 

accreted 

luminous 

satellites 

# 

surviving (<300kpc) 
luminous 

satellites (1O‘^M0) 

80% halo 

accretion 

time 

(Gyr) 

time of 

last >10% 

merger 

(Gyr) 

Number of 
particles 
in stellar 

halos 

HI 

115 

18 

4.16 

5.3 

8.3 

1908100 

H2 

102 

6 

3.26 

7.0 

9.2 

1890197 

H3 

106 

16 

4.17 

7.4 

8.9 

1689035 

H4 

97 

8 

4.11 

6.3 

8.3 

1628686 

H5 

160 

18 

2.40 

2.1 

10.8 

2531076 

H6 

169 

16 

2.45 

6.2 

10.5 

2873983 

H7 

102 

20 

2.44 

4.4 

7.4 

1797682 

H8 

213 

13 

2.58 

7.1 

9.3 

4246733 

H9 

182 

15 

2.52 

1.5 

10.0 

2528762 

HIO 

156 

13 

3.12 

2.9 

9.7 

2785418 

Hll 

153 

10 

2.56 

7.2 

9.0 

2598580 


TABLE 1 

Properties of the simulated stellar halos .Times refer to 

LOOKBACK TIMES. 


0.25 X 300 X 0.25 kpc^. Our models show that inner regions 
of the halos (r < 50 kpc) are consistently more homogeneous 
than the outer ones, both in terms of the spatial distribution 
of stars and in the overall chemical abundances. This is be¬ 
cause the inner halo builds up rapidly, with most of its stel¬ 
lar mass a lready in place in the first 4-5 Gyr of the Galaxy 
eyolution (iBullock & Johnstonl2005tlFont et alJ2005l) . and at 
these times the dynamical timescales were shorter and phase¬ 
mixing more efficient. The spati al distribution of the s t ars ha s 
been discussed in more detail by lBullock & .TohnstonI ll2005l) . 
who also note a similar behayior in the surface brightness and 
yelocity distribution of halo stars. Here we concentrate on the 
corresponding [Fe/H] and [a/Fe] distributions. 

We find that the inner ~ 50 kpc regions are, on ayerage, 
metal enriched ([Fe/H] > -1.2) and we trace this to the metal 
contribution of a few massiye (M* ~ satellites ac¬ 

creted early on and which contribute a significa nt fraction 
(40- 60%) to the total stellar mass of the halo (iFont et al.l 
l2005l) . The rapid growth of these satellites at early times fu¬ 
eled an intense star formation and therefore at the time of ac¬ 
cretion the satellites were already enriched in [Fe/H]. In the 
outer regions (r > 50 kpc) the streams are more spatially sep¬ 
arated and their yarious chemical enrichments are more dis¬ 
tinct. Both metal poor and metal rich streams are present, 
with the more metal poor ones originating in satellites with 
low star formation rates (these are typically low mass satel¬ 
lites). A more detail ed d iscussion of radial distributions of 
stars will be giyen in 33.31 

Figure |2l shows the ayerage [Fe/H] and [a/Fe] for all our 
Milky Way-type halos. The general trend in our simula¬ 
tions is that the more massiye a stellar halo, the more chem¬ 
ically eyolyed it is, i.e. the higher [Fe/H] and lower [a/Fe] it 
has (these relations are independent of the cutoff radius as¬ 
sumed for Mtfir)). The [Fe/H]-M* trend in Figure |2] is 
consistent with haying the same slope as that of the well- 
known relation between metallicity and luminosity seen in 

—2/5 

low-mass galaxies, [Fe/H] cx M* ' (dashed line). We note 
that our chemical feedback model is normalized such that 
suryiying satellites at z = 0 match the obsery ed metallicity- 
luminosity relation for Local Group gala xies (lLarsonlll974t 
iDekel & Siltd 1198^ iDekel & Wool 1200.31) . It is interesting 
howeyer that a similar relation propagates to the stellar halos 
themselyes. This follows from the fact that the slope of the 
relation for satellite galaxies is goyerned by their parent halo 



Fig. 1.— projections of chemical abundances in halos HI — H8 (in 

order from top left to bottom right). For each halo we show side by side 
the star particles color-coded in [Fe/H] (left sub-panels) and in [a/Fe] (right 
sub-panels). The maps span 300 kpc on the size. Each pixel in these maps 
extends 0.25 kpc in the x and z directions and 300 kpc in the y direction (i.e., 
in depth). The [Fe/H] and [a/Fe] values are mass weighted averages over all 
stars contained in each pixel. 

potential well depths. Therefore disrupted satellites also fol¬ 
low a [Fe/H] oc relation (see Figure 4 lFont et alJ2005ft . 
but with a lower [Fe/H] normalization because the se systems 
tend t o be accreted earlier a nd are less eyolyed dFont et alJ 
l2005t iRobertson et aljl2005l) . The same relation propagates 
to a stellar halo mass-metallicity relation because stellar halo 
mass correlates stro ngly with the number of disrupted massive 
satellites (see 33.2> . The more massiye disruptions a stellar 
halo contains, the more metal rich the stellar halo becomes. 
Similarly, because its progenitors are massiye, they are more 
metal rich, and the halo is higher metallicity as a result. 

Our results suggest that the hierarchical assembly is an 
important factor in determining the metallicity - mass rela¬ 
tion for stellar halos. A similar argument can be extended 
to qualitatiyely underst and the obseryed trend presented by 
iMouhcine et al.l d2005lJ) . who show that low luminosity spi¬ 
ral galaxies haye more metal poor central stellar halos than 
do high luminosity spirals. This qualitatiye trend can be 
explained in a hierarchical scenario: low-luminosity spirals 
tend to haye accretion histories dominated by correspondingly 
lower-luminosity dwarfs. Lower luminosity dwarfs are more 
metal poor and thus produce lower metallicity stellar halos. 
While one can deriye the [Fe/H] -M, relation based on simple 
scaling arguments that relate the metal loss to the depth of the 
centra l galaxy’s dark mater halo (e.g. iLarsonin 974t iTinslevI 
fl^ , this does not necessarily imply that the relation forms 
in situ. 

We caution that for a complete understanding of the ori¬ 
gin of the relation (both of its shape and normalization), one 
needs a more accurate treatment of gas physics processes in 
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Fig. 2.— Mass weighted [Fe/H] and [a/Fe] average values versus stellar 
mass for halos Fll — Fill. Averages are calculated within r < 50 kpc, a 
cutoff value which approximates the inner few tens of kpc currently probed 
by observations. The dashed line in the top panel corresponds to the [Fe/H] ~ 

fit of IDekel & Wod 1200311 . Errorbars represent weight averaged 25% 
and 75% values of the absolute spread in abundance ratios. 

satellites after they are accreted onto the main halo. These ef¬ 
fects are not currently included in our models, but we expect 
them to be of second order to the hierarchical mass assembly. 

3.2. Mass accretion histories and halo metallicities 

We can verify the hierarchical origin of the [Fe/H] -M* re¬ 
lation by investigating the mass accretion histories of the stel¬ 
lar halos. Given that in a typical mass accretion history the 
few massive satellites contribute significantly to the final stel¬ 
lar mass, we expect the overall halo metallicity (and hence the 
[Fe/H] -M* trend) to be determined mainly by these massive 
accretions. The time of accretion of massive satellites is also 
a related factor - satellites accreted later are likely to be more 
metal rich. 

Figure 0 plots for each stellar halo in our sample the 
lookback time of accretion of the most massive progenitors, 
specifically an average over those satellites which contribute 
more than 10% to the total stellar fraction. The general ten¬ 
dency is that more massive halos assemble more recently. Ad¬ 
ditionally, for halos of the same mass the spread in [Fe/H] can 
be also explained by differences in their formation time (e.g. 
the three halo models with M* ~ 1.5 x lO^M©). Therefore a 
protracted assembly of the higher mass stellar halos seems to 
be at the origin of the [Fe/H 1 -M* relation (see also similar 
results of iRenda et al.l20()5h . 

A Test Case: Milky Way and M31 

We further illustrate the relation between mass accretion 
histories and the overall metallicities with the case of Milky 
Way and M31. These are two large galaxies with roughly the 
same luminosity, but the Milky Way halo is more metal poor 
than that of M31, perhaps by as much as 1 dex. The metallic¬ 
ity distribution functio n (MDF) in the Milky Way halo peak s 
around [Fe/H] = -1.5 (iT.aird et aklll^ ICarnev et al.lfl9^ . 
whereas the MDFs measured at several locations in the inner 
“halo” of M31 peak at values between [Fe/H] -0 4 an d -0.7 
jReitzel & Guhathakurta|l^02t iBellazzini et al.ll2003l) . Al¬ 
though, as discussed below, the outer halo of M31 at [Fe/H] ~ 
-1.2 may be a more realistic estimate of the true underly¬ 


Fig. 3.— The lookback time of accretion of the most massive satellites, 
calculated as the average taccr of those satellites which contribute > 10% 
to the total mass of the stellar halo) versus stellar mass of the inner halo, 
M*(50kpc). Star symbols correspond to the most metal rich halos, triangles 
to halos of intermediate [Fe/H] and circles to those of lowest [Fe/H] in our 
sample. 

ing ste llar halo component of this galaxy llGuhathakurta et all 

Figure shows the mass accretion patterns, expressed as 
the fractional contribution of accreted satellites to the total 
stellar mass of the inner halo, for two halo models similar to 
inner Milky Way and M31: halo H5, with < [Fe/H] > ~ -1.3, 
and halo H4, with <[Fe/H]> ~ -0.9 (these are also the most 
extreme cases in our sample). The more metal poor halo has 
one major accretion of stellar mass, M* ~ accreted 

about 11 Gyr ago, while the more metal rich one has two 
major accretions, both of masses M* ~ IO^Mq and accreted 
about 8.5 Gyr ago. Thus the most massive, metal rich halo had 
an assembly history which includes more massive progenitors 
and which were accreted later. 

The inference from these results is that the difference in 
metallicities between the two galaxies can be explained if 
M31 experienced one or two more massive accretions than 
the Milky Way and/or a protracted merger history. 

However we note that the metallicity measurements in both 
Milky Way and M31 are not yet completely reliable. One rea¬ 
son is that most observations probe only the inner few tens of 
kpc and therefore may not be representative of the entire halos 
or may include additional galactic components. For example, 
recent observations in M31 extending further out into the halo 
find [Fe/H] ~ -1.2 for stars beyond 60 kpc, suggesting that the 
high metallicities measured in the inner regions include the 
contribution from M31 ’s e xtended bulge llGuhathakurta et all 
1200^ iKalirai et alJ 12001 Gilbert et al., in prep;). Addi¬ 
tional contamination may come from the recently discovered 
disk-l i ke component ex tending out to ~ 70 kpc llTbata et alJ 
1200^ iR^^hT^t^lTOOl There are also differences in the 
measurement methods. In the Milky Way metallicities are 
determined spectroscopically, while in M31 they are mostly 
determined p hotometrically (but see new spectroscopic m ea¬ 
surements of llbata et alJl2004 iGuhathakurta et al.ll2005bl) . It 
is therefore unclear whether the metallicity discrepancies be¬ 
tween the two galaxies are real, but ongoing wide-held spec¬ 
tral studies in these galaxies may be soon able to answer this 
problem. Nevertheless, our results are more general and sug- 
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Fig. 4.— Mass accretion histories expressed as the stellar fraction /* of 
the inner (r < 50 kpc) halo versus the stellar mass M* of the contributing 
satellites. The top panel corresponds to halo H5, a metal poor halo with 
<[Fe/H]> ~ -1.3, and the bottom panel to halo H4, a metal rich halo with 
<[Fe/H]> ~ —0.9. Empty circles show individual contributions of satellites 
to the stellai* fraction, and dashed lines their cumulative contribution. 

gest ways in which massive accretion events can increase the 
average metallicities of halos (eg. 1 massive dwarf galaxy, 
M* ~ IO^Mq, accretion is sufficient to increase the average 
metallicity of a Milky Way-type halo by about 0.5 dex). 

3.3. Radial Distributions 

As explained in the Introduction, metallicity gradients are 
a powerful tool for constraining the timescale of formation of 
galaxies. An often-cited hypothesis is that hierarchical forma¬ 
tion will result in a negligible metallicity gradient. In detail, 
however, some variations may occur, particularly in regions 
which continue to assemble until recent times. 

The radial distributions of [Fe/H] and [a/Fe] abundances 
for all simulated halos are shown Figures |3 and |6l The dis¬ 
tributions do not show significant global gradients, but many 
do show variations over scales of a few tens of kpc. 

The lack of large-scale gradients can be explained by ex¬ 
amining the mass accretion patterns of our simulated halos. 
The bulk of the stellar halos forms early, both in the inner 
and outer regions. The outer halos continue to accumulate 
stellar debris up to more recent times (~ 5 - 6 Gyr ago), but 
this represents only a small fraction of the total stellar mass 
dBullock & Johnstonl2005tlFont et alJ2005l) . Therefore, even 
if the later assembly of the outer halo favors an influx of chem¬ 
ically enriched material, this is not sufficient to create sig¬ 
nificant large-scale gradients. Moreover, the chemical abun¬ 
dances are directly correlated with the ages of the stars, rather 
than with the time of accretion. Figure0plots the distribution 
of lookback accretion times (left panels) and the ages of stel¬ 
lar populations (right panels) versus the radial distance, for 
all halo models. As this figure shows, even if the stellar mate¬ 
rial is accreted recently, it usually contains a mixture of both 
young and old stars that dilutes the metallicity differences. 

Many of our halos do show distinct small-scale abundance 
gradients. For example, in halos H5, H6 and H9 the aver¬ 
age [Fe/H] increases with r in the inner 100 kpc (similarly, 
[a/Fe] decreases over the same distance). The gradients in 
these halos arise from the spatial distribution of their stellar 
populations. Figure 0 shows that the inner regions of these 


three halos contain mostly old (^ 10 Gyr) stellar populations, 
whereas the outer regions have on average stellar populations 
a few Gyr younger. This result has implications for searches 
for the lowest metallicity stars in the stellar halo and compli¬ 
mentary attempts to use the metallicity distribution function 
to constr ain the nature of reionization and PopIII s tar forma¬ 
tion (e.g. lTumlinson et alJl2004HBeers et alJl2005l) . Specifi¬ 
cally, the lowest metallicity stars are expected to inhabit the 
most central regions of the Galaxy (note the inner-most ra¬ 
dial bin in Figure |^. Thus general attempts to use observed 
metallicity distribution functions to constrain early chemical 
enrichment must take special care to account for spatial biases 
in observational samples of stars. 

The radial distributions of accretion times, ages and chemi¬ 
cal abundances give an insight into the composition of stellar 
populations in the simulated halos. Our results suggest that 
stellar halos of galaxies like the Milky Way should contain 
a large fraction of old (> 10 Gyr) stellar populations which 
were accreted in the first few Gyr of the mass assembly. In 
addition to the underlying population of old stars, some ha¬ 
los contain small fractions of intermediate age (~ 7-10 Gyr) 
populations. Irrespective of their age, stars can be either metal 
poor or metal rich, depending of the mass and star formation 
rate of their progenitor satellite. 

Interestingly, in terms of its age and metallicity distribu¬ 
tions, the Milky Way halo lies at the edge of our sample sim¬ 
ulated halos. The stars in the Milky Way halo are mostly old 
and metal poor, whereas in our models a significant number 
of old stars are already metal enriched. (Note however that 
the typical observations cover only the inner ^ 10-20 kpc of 
the halo and in this range some of our halos are old and metal 
poor). Our results seem to suggest that the Milky Way halo 
is an atypical case in the range of galaxies of similar mass. 
The metallicity constraints rule out an accretion of a massive 
(M* ~ IO^Mq) halo progenitor at early times when the bulk of 
the halo formed. We also note that recent observations show 
that many other bright spiral galaxies of similar mass have 
more metal rich halos than that of the Milky Way. A more 
detai led d iscussion of these observations will be given in Sec¬ 
tion 34.11 

3.4. Mining the Halo for Stellar Streams 

Tidal streams are fossil records of past accretions and there¬ 
fore can be used to reconstruct the merger histories of galax¬ 
ies. Their detection however is difficult due to their low sur¬ 
face brightness and to the loss in spatial coherence the longer 
they orbit in the Galaxy. The identification of stellar streams 
can be improved by extending the parametric space to include 
new dimensions like the kinematic parameters and chemical 
abundances. 

3.4.1. Kinematics 

If kinematic data is available from observations, one can use 
the (£, Lj, L) space to identify stellar streams. This paramet¬ 
ric space is particularly useful for epochs of slow evolution, 
when the time required for stars to exchange energy and angu¬ 
lar momentum is long compared with the age of the Galaxy, 
and E, L, and L can be approximated as invm ants of mo¬ 
tion (e.g. fEggen, Lynden-Bell & Sandagell962h . In this case, 
the tidal debris from different accretions is expected to clump 
in distinct locations in this space. Numerical simulations of 
satellite accre tions in a fixed potential c onfirm that this is in¬ 
deed the case (iHelmi & de Zeeuwl200r)h . 
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Fig. 5.— Radial distributions of [Fe/H] (left panels) and [a/Fe] (right panels) for halos HI — H5. The chemical abundances are weight averages in radial shells 
of width dr = 10 kpc. En'or bars represent weight averaged 25% and 75% values of the absolute spread in abundance ratios within each radial bin. 


In contrast with previous studies we model the mass ac¬ 
cretion history in a cosmological context and include a time- 
evolving gravita tional potential which grow s with the mass of 
the Galaxy (see iBullock & .lohnstonll2005[ for details). We 
want to see if, with these additions, the tidal debris can still 
be recovered in the {E,L^,L) space. Figure|^plots the dis¬ 
tribution of halo HI stars in the L--L, E - L~ and E-L 
planes at different ranges of lookback accretion times of the 
progenitor satellites. This shows that the recent accretions 
(taccr <7-9 Gyr) occupy relatively distinct locations in these 
planes and their identification is therefore possible. 

The clumps in the low E and L regime are more difficult 
to separate, in part because of their superposition, but also 
because this is also the regime of rapid growth of the galaxy 
when the energy and angular momentum of stars change sig¬ 
nificantly. Additional effects like the interaction between 
satellites and the non-axisymmetric oscillations in the Milky 
Way potential, not mode led in our study, w ill further disperse 
the stars in phase-space llMaver et aU2002h . 

The kinematic identification of stellar streams in the halo 


requires high accuracy measurements of proper motions over 
extensive areas of the sky. Future astrometric Galactic surveys 
such as GAIA will be able to provide these data and identify 
some of these clumps llHelmi & de Zeeuwl200flt) . 

3.4.2. Chemical Abundances 

Chemical abundances can be used as an additional dimen¬ 
sion to the ( £, Ly, L)^^am^~ic space to in^sdg^e the origin 

rating the available kinematic parameters in different metal- 
licity ranges, some observational studies have found stars 
which stand out from the rest of the halo, implying a merger 
origin (e.g. ICarnev et alJ 11995 Maiewski, Munn & Hawle^ 
ri9^ . However, given that even stars from the smooth 
halo may originate in mergers, it is instructive to investigate 
whether chemical abundances can provide any new informa¬ 
tion about the merger history of the galaxy. 

In particular, chemical abundances may be able to constrain 
the merger timeline of the Galaxy. In Figure |3 we have used 
the time of accretion to separate the clumps in the {E,L,, L) 
space. This parameter is not easily measurable, but fortu- 








































Phase-Space Distributions of Chemical Abundances 


7 




<D 


<D 

ti-i 




q; 

tin 


p: 

0) 

Ijl- 



0) 

1^ 


a; 


<3 


0) 


pCH 


0) 

Plh 



H6 


H7 


H8 


H9 


H10 


H11 


r [kpc] 


50 100 150 200 

r [kpc] 


Fig. 6.— Radial distributions of [Fe/H] (left panels) and [oi/Fe] (right panels) for halos H6 - HI 1. The abundance ratios are weight averaged values within 
radial bins of 10 kpc. EiTor bars represent weight averaged 25% and 75% values of the absolute spread in accretion times and ages within each radial bin. 


nately, we can use the natural clock provided by the evolu¬ 
tion of Fe and a-elements. Each satellite galaxy is expected 
to have a distinct [a/Fe] - [Fe/H] evolution, however it has 
some generic features such as a plateau at high [a/Fe] and low 
[Fe/H] and then a steady decrease towards high [Fe/H] (both 
the length of the plateau and the slope of the decrease be¬ 
ing modulated by the star formation rate, hence the mass 
of the satellite). The plateau is maintained by the enrich¬ 
ment in a-elements in Type II supernovae operating on short 
timescales (of a few Myr) and the [a/Fe] decrease is driven 
by the later enrichment in Fe produced in Type la supernovae. 
Prolonged episodes of star formation will tend to decrease the 
[a/Fe] ratios of stars in a satellite. This suggests an associa¬ 
tion between the time of accretion of satellites and their over¬ 
all [a/Fe]. 

Figure |3 plots the weight averaged [a/Fe] versus 
[Fe/H] values for satellites accreted onto halo HI, high¬ 
lighting the times of accretion. This shows that the lowest 
[a/Fe] values are associated with recently accreted satellites, 
these systems having more time available to form stellar pop¬ 


ulations with [a/Fe] below solar. This is more clearly seen in 
Figurenol where we plot the distribution of stars in the E—L, 
and E-L planes, separated in different [a/Fe] ranges. Most 
of the stars with [a/Fe] < 0.05 are in the high energy and 
angular momentum regime, associated with recent accretions. 

This also suggests that by selecting an [a/Fe] < 0.05 cut, 
one can effectively remove the bulk of the old, well mixed 
halo and identify only the recently accreted streams. Fig¬ 
ure ^ illustrates this with a comparison between the surface 
brightness maps of some of the simulated halos, first by in¬ 
cluding all stars and alongside including only those stars with 
[a/Fe] < 0.05. With this cut, the central regions are effec¬ 
tively devoid of the well mixed streams and contain only the 
cold streams. 

How feasible is our selection method for detecting streams 
in the halo? Figure ^ shows that even without the [a/Fe] cut 
the tidal streams are faint. Only a few of these streams, typ¬ 
ically 1 or 2 per galaxy, can be detected with current capa¬ 
bilities which have a limiting surface brightness in V-band of 
/iy ~ 30-31 mag/arcsec^. The proposed method can be used 























8 


Font et al. 


10 

5 

F ■ - - • - M . I f * s f ^ J = i 

10 


5 


10 

5 

■ 1 1 1 1 1 1 l*Ai 1 1*1 1 1 1 1 1 1 

10 

ii.... . .... 

5 

^ . 11 i H J 

k , 1 , 1 , , , , 1 , , , , 1 , 1 , 

10 

5 

r iH.H ..IH 

10 

-Illilll i i 


i I I 


... , 1 . 

10 

5 

. 

= 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

10 

5 

, JlJ 

1 1 ■ * *|" 

10 

5 

= , , , , 1 , , , , 1 , , , , 1 , , , ,S 

10 


,—, 5 

^.* , , 

o 10 



- ■* '~ lllllllil ] 

O 5 



0 50 100 150 200 

r [kpc] 



[kpc] 


Fig. 7.— Radial distributions of lookback accretion times (left panels) and ages of stellai' populations (right panels) for halos HI — HI 1. The abundance ratios 
are weighted averages within radial bins of 10 kpc. Errorbars represent weight averaged 25% and 75% values of the absolute spread in abundance ratios in each 
radial bin. 


to highlight a few cold streams of about this limiting magni¬ 
tude in the inner halo. 

The prospects for detecting giant stars are also reasonable. 
Assuming the luminosit y of halo stars is similar to the o ne 
of globular cluster Ml2 iHargis. Sandauist & Roltell2004i . a 
IO'^Lq halo will contain a few 10® giants (out of these, 10® gi¬ 
ants will be 1 V-band mag above the horizontal branch (HB) 
or about 10"^ giants 2 V-band mag above the HB, respec¬ 
tively). According to our models, 1 - 10% of the stars meet 
the [a/Fe] < 0.05 cut. This implies, for example, that a halo 
survey of 10® giants with Imag above the HB will contain 
~ 10®- 10"^ giants in cold streams. 

4. COMPARISONS WITH OTHER STUDIES 
4.1. Stellar Halo Metallicities 

A series of recent observational studies have found evi¬ 
dence in support of the idea that the metallicity/age of the 
Milky Way stellar halo is peculiar. It seems that most 
other bright spir al galaxies are more m etal rich than the 
Milky Way (e.g. iMouhcine et"^ 12005^ . M31 is a prime 


example, but recent observations show that this is also the 
case for several spiral galaxies outside the Local Group 
ilMouhcine et al.l2005hl) . as it is for about 1000 spiral galaxies 
in t he SPSS sample, as inferred from th e average of their col¬ 
ors iZibetti, White & Brinkman~n1l2004li . Similarly, in terms 
of their stellar composition, the observations suggest that all 
halos have an underlying, uniform population of old metal 
poor stars and that most luminous galaxies have additional 
metal rich popula tions, which can b e either old or of inter¬ 
mediate age (e.g. lHa nis_etaLl [19991 |Harris_&JIarrislEQQ2l; 
Bellazzini et alJ^03l IMouhcine et alJl2005bt iGallazzi et alJ 

2nn.5h . As a bright galaxy but with an old and metal poor 
stellar halo, the Milky Way may indeed be somewhat rare. 
Indeed, none of our 11 simulated stellar halos have metallic¬ 
ities as low as that of the Milky Way <[Fe/H]>~ -1.5, 
although our most metal poor halo is close <[Fe/H]> ~ -1.3 
(see Figure |2}. Our interpretation would be that the Milky 
Way is unusually lacking in massive satellite mergers at early 
times. While rare, this would not be unthinkable within the 
expectations of ACDM since the average number of accreted 
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Fig. 8.— Distribution of stars in halo HI in the L space (left panels), 
in £ —L;; (intermediate panels) and in E-L (right panels). Different panels 
from top to bottom correspond to diffferent lookback accretion times of the 
satellites where stars originate from. 


Fig. 10.— Distribution of stars in halo H2 in the E—L- space (left panels) 
and in E—L (right panels). Top sub-panels show all stars, while the remaining 
three sub-panels from top to bottom correspond to different [oi/Fe] ranges. 
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Fig. 9.— [a/Fe] versus [Fe/H] for all baryonic satellites accreted onto halo 
HI. Top: Satellites are plotted with different symbols function of their look- 
back time of accretion: squares denote satellites accreted at taccr > 9 Gyr, 
triangles those accreted 5 < taccr < 9 Gyr; and star symbols those accreted 
taccr < 5 Gyr. Bottom: Small and large symbols denote whether satellites 
are fully disrupted or still containing some bound material at present time 
{t = 0), respectively. En'orbars in both represent 10% and 90% of the weight 
averaged values. 


massive satellites is small (IZentner & Bullockl2003l) . 

As shown in Figure |2j the metallicities of our simulated 
halos have a spread of les s than 1 dex. Th is disagrees 
somewhat with the results of iRenda et al.l il200.^ . who report 
that their simulated halos have a larger [Fe/H] spread, about 
1.5 dex, when restricted to have Milky Way-ty pe luminosi¬ 
ties. More specifically, the iRenda et al.l ll200-^ simulations 
produce Milky Way-type halos as metal poor as <[Fe/H]> ~ 
-2. We suspect that the difference arises from the fact that 



Fig. 1.— V-band surface brightness maps for halos HI, H2, H3, H4, H5, 
H7, H8 and HIO (in order from top, left to right). Left sub-panels show the 
surface brightness map of all stars and right sub-panels the corresponding 
surface brightness of stars with [a/Fe] < 0.05. All maps span 300 kpc on the 
size. 










































10 


Font et al. 


iRenda et alJ (l2005l) allow stellar halo stars to originate via in 
situ star formation in their simulations, while our stellar ha¬ 
los form entirely from accreted dwarfs. Thus the metallicity 
spread in observed halos may provide an interesting avenue 
for testing the idea of accreted vs. in situ stellar halo forma- 
tion. The smal l spread in stellar halo metallicities seen by 
iMouhcine et alJ ( l2005bl) may favor an accretion-dominated 
model, although the evidence is weak as of now. A sec- 
o nd difference may a rise from the cosmological conditions 
in iRenda et alJ ( l2005l) . which are set via a semi-cosmological 
top-hat sphere with standard CDM fluctuations, while ours 
rely on Extended Press Schechter within an ACDM frame¬ 
work. It is interesting to note that our restriction on accretion 
histories to be those without recent major mergers would (if 
anything) bias our distributions to more metal poor halos — 
just the opposite bias that w ould be required to match recon¬ 
cile our results with those of iRenda et al.l(l2005l) . 


4.2. Spatial Distributions 

Our models predict that tidal streams from massive satel¬ 
lites should stand out in stellar overdensities and be more 
metal rich. The available observational evidence seems to 
support our finding. A wide survey of the inner M31 halo 
reveals an increase in metallicity in fields associated with the 
giant stellar stream JFerguson et alJ2d(^ . Numerical model¬ 
ing of the stream finds t hat the progenitor satellite was indee d 
massive, M* ^ Mq dFont et al. 1120041 fFardal et al.ll2n0-5h . 
Inhomogeneities in color magnitude diagrams are also seen in 
the M31 h alo and are believed to be associated wit h satellite 
accretions liRrown et al.l2Q03irFerguson et al.l2005l) . 

While metallicity inhomogeneities are commonly detected, 
the observational evidence regarding radial gr adients is so far 
inconclusive. Stu dies of globular clusters (e.g. lSearle & Zinr3 
WM inn|122^or of low metallicity stars in the Galaxy 
jBeers et alJboOSt) do not show obvious gradients in [Fe/H]. 
However, some studies report a decline of \ ct/Fel ratios of halo 
stars with radial distance in the Galaxy ilNissen & .Schiisteii 
Il997h . In M31, measurements in several fields extending up 
to 40 kpc in the halo are co nsistent with a small, less than 
1 dex, [Fe/H] radial gradient JBellazzini et al.l2003 ). 

Our results suggest that metallicity gradients of less than 
1 dex over distances of a few tens of kpc are consistent 
with the hierarchical structure formation scenario. How¬ 
ever, for a more accurate comparison between our models 
and observations, we have to take into account the meth¬ 
ods of measurement used in observations. For example, in 
external galaxies radial distances cannot be currently deter¬ 
mined. What observations generally measure are the separa¬ 
tions from the centers of galaxies, which are projected radii 
on the plane of the sky. Moreover, current observations probe 
smal l fields in the stellar halos at different projected radii 
(e.g. iBellazzini et^l2003ft . We illustrate this with a mock 
observation of one of our halos (HI). Figure [El shows the 
[Fe/H] and [a/Fe] distributions versus the projected radius R 
on the plane (x,z). The large symbols are averages in cylindri¬ 
cal shells, while the small symbols show the decomposition of 
these values at equidistant azimuthal angles in each cylindri¬ 
cal shell. Thus each small symbol approximates a projected 
field of observation a few kpc on the size, similar to the size 
of observational fields. Figure El demonstrates that a large 
spread in either [Fe/H] or [a/Fe] may exist at a given R and 
therefore a random sampling with small fields over different 
projected radii can result in a detection of spurious abundance 
gradients. Our results therefore caution that a wider spatial 




Fig. 12.— [Fe/H] and [o/Fe] distributions versus cylindrical radius R in 
Halo HI, (R = V-F + z^). Large pentagons correspond to weight averaged 
[Fe/H] or [cr/Fe] values in cylindrical shells of width dR = 10 kpc. Small 
squares correspond to weight averaged [Fe/H] or [cr/Fe] values at equidistant 
azimuthal angles </< £ [0, Ivr] in a given cylindrical shell. (Note that the num¬ 
ber of angles </< probed in each shell increases with R, so as to preserve the 
size of the mock observational fields). 

sampling of halos is necessary in order to confirm the current 
detections of abundance gradients. 

5. CONCLUSIONS 

The motivation of this work was to provide predictions of 
the hierarchical formation paradigm that can be tested with 
the data from current and upcoming astrometric and spectro¬ 
scopic observations. In this study we have explored the global 
properties as the phase-space distribution of the [Fe/H] and 
[a/Fe] chemical abundances in a series of Milky Way-type 
stellar halos formed in a ACDM Universe. Our results can 
be summarized as follows: 

- For a fixed dark matter halo mass, we predict that stellar 
halo metallicities increase with stellar halo mass. The resul¬ 
tant [Fe/H] -M^ relation for stellar halos arises because high- 
mass stellar halos tend to have accreted a larger number of 
massive satellite progenitors. The relation propagates hierar¬ 
chically because massive progenitors are more metal rich than 
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their low-mass counterparts. 

- Chemical abundance distributions provide important links 
to the minor merger histories of galaxies. For example, the ob¬ 
served metallicity differences between the (metal poor) Milky 
Way and (metal rich) M31 stellar halos may arise because 
M31 destroyed one or two more massive satellites than the 
Milky Way. 

- In agreement with observations, we find that all simu¬ 
lated halos contain an underlying population of old metal poor 
stars. More massive halos contain additional intermediate age 
populations accreted more recently. 

- In contrast with the Milky Way, some of our halo mod¬ 
els contain old stellar populations which are metal rich, origi¬ 
nating in massive satellites (M* ~ 10®Mq) accreted early on. 
The metallicity constraints seems to suggest that the Milky 
Way has an atypical mass accretion history, lacking in mas¬ 
sive early accretion events, compared with other galaxies of a 
similar mass. 

- We predict metallicity gradients in [Fe/H] and [a/Fe] in 
our hierarchically-formed stellar halos that are often non- 
negligible. Gradients extend typically over a few tens of kpc 
and may be as large as 0.5 dex in [Fe/H] and 0.2 dex in [a/Fe]. 

- We predict that most metal poor stars in the Galactic halo 
are buried within the central ^ 5kpc of the Galaxy. This will 
likely be important for efforts to use the observed count of low 
metallicity stars to constrain models of cosmic reionization 
and the transition from Pop III to Pop II star formation (e.g. 
iTumlinson et alJ2004) . 

- We suggested a method to identify cold tidal streams 
in stellar halos by selecting stars with [a/Fe] < 0.05, which 
trace streams from recently accreted satellites. If used to 
complement the common kinematical identification methods, 
the proposed selection criterion may improve the detection of 
tidal streams. 

Combining kinematic and chemical abundance information 
proves to be a powerful tool for understanding how galaxies 
form and evolve. In anticipation to the wealth of data to be 


Abadi, M. G., Navarro, J. R, & Steinmetz, M. submitted to MNRAS 
I astro-ph/0506659 i 

Altmann, M., Catelan, M. & Zoccali, M. 2005 A&A, 439, L5 
Beers, T. C. et al., to appear in “From Lithium to Uranium: Elemental Tracers 
of Early Cosmic Evolution", lAU Symposium 228, V. Hill, P. Erancois & 
E. Primas, eds. i astro-ph/0508423 i 
Bekki, K. & Chiba, M. 2001, ApJ, 558, 666 

Bellazzini, M., Cacciari, C., Federici, L., Fusi Pecci, F, Rich, M. 2003 A&A, 
405, 867 

Bland-Hawthorn, J. & Freeman, K. 2000, Science, 287, 79 
Blumenthal, G. R., Faber, S. M., Primack, J. R. & Rees, M. J. 1984, Nature, 
311, 517 

Brodie, J. & Huchra, J. P. 1991, ApJ, 379, 157 

Brook, C. B., Kawata, D., Gibson, B. K., & Flynn, C. 2003, ApJ, 585, L125 
Brook, C. B., Kawata, D., Gibson, B. K., & Flynn, C. 2004, ApJ, 349, 52 
Brown, T. M et al. 2003, ApJ, 592, 17L 

Bullock, J. S. & Johnston, K. V., ApJ, accepted (astro-ph/0506467 i 

Bullock, J. S., Kravtsov, A. V., & Weinberg, D. H. 2001, ApJ, 548, 33 

Carney, B. W., Laird, J. B., Latham, D. W., Aguilar, L. A. 1996, AJ, 112, 668 

Chiba, M. & Beers, T. C. 2000, AJ, 119, 2843 

Dekel, A. & Silk, J. 1986, ApJ, 303, 39 

Dekel, A. & Woo, J. 2003, MNRAS, 344, 1131 

Diemand, J., Madau, R, & Moore, B. 2005, MNRAS, in press 

Dinescu, D. 1. 2002, ASPC, 265, 365 

Eggen, O. J., Lynden-Bell, D. & Sandage, A. R. 1962, ApJ, 136, 748 
Eardal, M. A., Babul, A., Geehan, J. J., & Guhathakurta, R, submitted to 
MNRAS (astro-ph/0501241 I 


delivered by current and future Galactic surveys, it is impor¬ 
tant to develop and improve on numerical simulations of the 
chemical evolution of the Galaxy in a realistic cosmological 
context - such as those presented in our study - and to devise 
criteria to identify more tidal streams. 

While our methods successfully reproduce the gross chem¬ 
ical abundance properties of the Galactic halo and satellite 
dwarf galaxies, our treatments of dissipational physics includ¬ 
ing star formation, gas cooling, and energetic feedback mech¬ 
anisms may be too simplistic to recover the detailed proper¬ 
ties of the modeled systems with high accuracy. For instance, 
the assumed truncation of star formation in dwarf systems af¬ 
ter they accrete into their host galaxy is very approximate. 
Gravitational torques applied by the central galaxy and other 
satellite systems may induce episodes of efficient star forma¬ 
tion in the dwarfs we model, especially near pericentric pas¬ 
sage. Such bursts occurred in the LMC and SMC, as inferred 
from their stellar popula tions (e.g. lSmecker-Hane et al.l2002t 
iHarris & 7aritskvll2004) . and should be modeled if more ac¬ 
curate chemical enrichment histories for dwarf systems are to 
be achieved. 

These complications speak to the need for high-resolution 
hydrodynamical simulations of the formation of the Galac¬ 
tic stellar halo and dwarf satellite population. While recent 
hydr odynamical simula tions have begun to address this issue 
(e.g. lBrook et al.ll2004ll . they have not yet matched the mass 
resolution necessary to capture the entire population of dwarfs 
galaxies we follow in our collisionless simulations nor the full 
complexity of our chemical enrichment model. We acknowl¬ 
edge the importance of the previous work and look forward to 
increasingly sophisticated hydrodynamical simulations of the 
cosmological formation of the Galactic stellar halo. 
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